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Abstract 

Multiple-breath-washout (MBW) measurements are regarded as a sensitive 
technique which can reflect the ventilation inhomogeneity of respiratory air- 
ways. Typically nitrogen is used as the tracer gas and is washed out by 
pure oxygen in multi-breath-nitrogen (MBNW) washout tests. In this work, 
instead of using nitrogen, ^He is used as the tracer gas and a multiple-helium- 
breath-washout (MBHW) system has been developed for the lung function 
study. A commercial quartz tuning fork with a resonance frequency of 32768 
Hz has been used for detecting the change of the respiratory gas density. 
The resonance frequency of the tuning fork decreases linearly with increas- 
ing density of the surrounding gas. Knowing the CO2 concentration from 
the infrared carbon dioxide detector, the helium concentration can be de- 
termined. Results from 12 volunteers (3 mild asthmatics, 2 smokers, 1 with 
asthma history, 1 with COPD history, 5 normal) have shown that mild asth- 
matics have higher ventilation inhomogeneity in either conducting or acinar 
airways (or both). A feature has been found in single breath washout curve 
from 4 smokers with different length of smoking history which may indicate 
the early stage of respiratory ventilation inhomogeneity. 
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1. Introduction 

1.1. Physiological Background 

Multiple-breath washout (MBW) is regarded as a sensitive technique to 
study the ventilation inhomogeneity in conducting and acinar airways [? ? 
?????]. In the majority of MBW measurements, subjects breathe 
in pure oxygen gas, replacing the air and washing nitrogen gas (the tracer 
gas) out of lungs breath by breath. The concentration of the tracer gas 
along with the respiratory tidal volume is carefully monitored as shown in 
figure 1. By plotting out the concentration of the tracer gas versus the 
expired tidal volumes, each single breath can be divided into three phases as 
shown in the inset of figure 1 uppergraph. Phase I corresponds to the pure 
inspiratory gas and there is no tracer gas. Phase II, with a sharp increase 
of gas concentration, is the transition between the inspiration and the gas 
remaining in the lung [? ]. Phase III corresponds to the alveolar plateau 
with a positive slope [? , pp.296]. 

The analysis of MBW focusses on phase III. The mechanism which causes 
the phase III slope to be slightly positive is inhomogeneous ventilation due to 
the asymmetrical structure of alveolar airways and the uneven distribution 
of alveoli. The inspired gas with lower concentration is mixed with the gas 
remaining in alveolar sacs and is distributed unevenly to different units of 
the lungs [?????]. After an inspiration, the gas mixture in those 
respiratory units with less ventilation, (i.e., more difficult for the inspired 
gas to diffuse into,) has a higher concentration of tracer gas and is flushed 
out less easily. The gas mixture in the better- ventilated units relatively has a 
lower concentration of tracer gas since the gas is diluted well and contributes 
to the early part of the alveolar expiration. Less- ventilated regions contribute 
to the later part of the phase III which causes the positive phase III slope. 
The washout results from individuals with lung diseases have steeper phase 
III slopes compared to normal subjects because of the greater inhomogeneity 
of their respiratory ventilation [?????]. Other minor effects such as 
body posture or inspiration and expiration flow rates can also increase the 
phase III slopes [? ? ]■ 

1.2. Physiological Analysis 

1.2.1. Lung Clearance and Two Compartment Model 

During washout, for unevenly ventilated lungs, the alveolar gas concen- 
tration becomes unequal from unit to unit and the flrst breath in the washout 
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Figure 1: The typical nitrogen washout curve measured in University of Leicester from 
a volunteer with a commercial nitrogen washout system (Ultima PF, Medical Graphics 
Cooperation, USA) . The upper graph shows the concentration of nitrogen and the lower 
graph shows the respiratory tidal volumes during the washout process. The arrow indicates 
the start of inspiring pure oxygen. The inset in the uppergraph shows a single breath curve 
where the nitrogen concentration is plotted against the tidal volume. I, II and III represent 
phases I, II and III. S7/7 represents the phase III slope which is fitted by a straight dashed 
line. 
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should reflect the gas flow the most from the better-ventilated units. The 
degree of ventilated inhomogeneity arises progressively less from those better- 
ventilated units and more from the poorly- ventilated units in the second and 
subsequent breaths. This is because the alveolar gas is washed out more 
quickly from the better-ventilated units than from poorly-ventilated units. 

For an ideal container with a perfect ventilation, i.e., the inspired dilut- 
ing gas mixes perfectly well with the tracer gas remaining in it, the mean 
concentration of expiration is exponentially related to the breath number in 
the MBW test with constant tidal volumes. The relation can be given by [? 



where n is the breath number, c„ and cq are the tracer gas concentration 
from the nth breath and from the lungs before washout starts, FRC is the 
residual capacity (or the volume of the container before inspiration), and Vt 
is the volume of tidal breath. Equation 1 can be written as 



Hence for the ideal lung with steady respiratory tidal volumes, the log- 
arithm of the mean concentration is linearly related to the breath number. 
However, for real lungs, the uneven distribution of inspired gas is washed out 
asynchronously and it results in a nonlinear relationship between logCn and 
n. The interpretation of the washout curve from MBNW was first studied 
by Robertson et al. [? , pp. 312] [? ]. The washout curve corresponds to the 
sum of two (or more) exponential functions representing two (or more) com- 
partments independently ventilated at different rates as shown in figure 2(b). 
For patients with airway diseases, the respiratory ventilation inhomogeneity 
reflects in a more curved hne [? ] . 

The lung contains units with different ventilation rates and the tidal vol- 
ume also varies breath by breath. In this study, the lung is simply mod- 
elled as two compartments with different ventilation rates leading to lung 
clearance curves that are fitted by the sum of two exponential curves with 
respect to the turnover TO [? ]. The turnover is defined as the cumu- 
lative expired tidal volumes divided by functional residual capacity FRC, 
i.e., TOn = YIi^t/FRC. The decay rates and the volume ratio of the two 
compartments are also compared for different subjects. Equation 1 is then 
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Figure 2: (a) The perfect lung model as a well-ventilated container with perfect gas 
mixing. The multiple-breath nitrogen washout curve is a straight line, (b) Two (or more) 
compartment model with two different ventilated rates. The nitrogen washout curve is a 
sum of two exponential curves [?][?, pp.312]. 
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Considering the two compartments 1 and 2 with the capacities FRCi, 
FRC2 and tidal volumes T^, V2, the lung clearance mechanism can be given 

by 

FRC = FRC1 + FRC2, (4) 

Vt = V^ + V2, (5) 
n I o ^ FRCi FRC2 . s 

'^' + ''^^ = FRC + FRC^'^'^^FRC^ + FRC^''"'' 

( FRCi \™ jjTO (rr\ 
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where Cn=o,i,2, - n-i are the mean concentration of gas and VT,n(n=o,i,2, - n-i) 
are the tidal volumes from the initial gas in the lungs, the first, second, and 
{n — l)th expiration, respectively. Ri and R2 represent the washout rates 
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from those two compartments. The more details of the calculation is given 
in Appendix A. 

The lung clearance index {LCI) has been commonly used as a measure 
of lung health and is defined as the total expired volume during washout 
divided by the FRC, for a fall in tracer gas concentration of a factor 40, i.e., 
the maximum turnover of the gas concentration above about 2% [???]. It 
has been shown that the LCI is more sensitive than the spirometry in cystic 
fibrosis (CF) adults ,and its mean value has been measured 13.1±3.8 in CF 
patients and a much lower mean value of 6.7±0.4 in normal subjects [? ? 
??????]. Other studies have also shown the increasing LCI with 
the age [? ] and a higher range of LCI in emphysematous patients [? 
]. In this study, we compare the LCI values from 12 volunteers; however, 
the LCI values have not been found to differ significantly between the mild 
asthmatics or smokers and the normal subjects. 

1.2.2. Two Indices Scond and Sadn 

Each breath of the MBW can be regarded as a single washout curve 
composed of phase I to III. By analysing the phase III slopes (Sjjj) from 
each breath, the washout results provide two indices Scond and Sadn which 
reflect the degree of ventilation inhomogeneities in conducting and acinar 
airways, respectively. For a perfect ventilated lung, the phase III is flat 
and Sjji is near zero. In real lungs, Sju is not only positive but also the 
normalised value, i.e., Sju devidcd by the mean tracer gas concentration, 
increases breath by breath. Scond value is taken as the increasing rate of the 
noramhsed phase III slopes after the flrst few breaths. The higher Scond value 
reflects the greater ventilation inhomogeneity in conductiing airways [? ? ? 
? ? ]. Sacin is the nomralised phase III slope from the flrst breath with the 
subtraction of Scond, i-G-) 

Sacin = Stii{1) — Scond ' TO (9) 

The higher Sadn value the greater ventilation inhomogeneity in acinar air- 
ways. More details of the physiological theory on Scond and Sadn can be 
found in [???]. 

Figure 3 shows the normallised phase III curves of helium washout mea- 
surements from a normal and a mild asthmatic volunteer. The normalised 
phase III slopes {Sm) of the asthmatic subject is apparently higher than of 
the normal subject as well as the increasing rate of the Sm, i.e., the larger 
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Figure 3: The helium washout resuhs of normahsed phase III slopes versus turnover from 
a normal subject and a mild asthmatic volunteer. 

Scond value, which imphes the higher ventilation inhomogeneity in the con- 
ductive airways. 

1 . 3. Multiple- Breath- Washout Measurements 

Nitrogen is usually chosen as the tracer (MBNW) because it is the gas we 
normally breathe and has no direct influence on physiology unlike oxygen. 
SFg is also becoming more frequently used in low concentrations [? ? ? ] but 
in this study, ^He gas is used as the tracer gas (MBHW) and is washed out by 
room air. Since helium not only has smaller density which is frequently mixed 
with oxygen to aid breathing but also has higher diffusivity than nitrogen, it 
is belived to be able to reach deeper into our lungs in a given time. Therefore, 
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helium washout may be a more sensitive maneuver reflecting the ventilation 
inhomogeneity especially in small airways compared to nitrogen washout. 
The other important point is that '^He docs not diffuse through the airway 
wall into the blood vessels but is washed out of the lungs by air. 

In our MBHW system, a quartz tuning fork (QTF) with a resonance fre- 
quency 32768 Hz is used as the gas density sensor. The resonance frequency 
of the QTF is linearly related to the surrounding gas density [???]. 
The helium concentration is calculated from the tuning fork signal having 
subtracted the contributions to the density of all other gas components, in- 
cluding carbon dioxide, oxygen, nitrogen, and water vapour. In our system, 
the carbon dioxide is detected by the infrared sensor and the water is flltered 
out by an ice bath. The oxygen concentration is given by the assumption of a 
constant gas exchange rate from the pre-washout measurement. The helium 
washout curve is shown in figure 4. 

We present results on 14 volunteers to demonstrate typical signals and to 
illustrate the information that can be derived from the MBW. We have per- 
formed the washout measurements on three mild asthmatics, four smokers, 
two having mild asthma or COPD history, and six five subjects. The peak 
expiratory fiow is also measured for each subject. 

2. Materials and Methods 

2.1. Principle of Measurement 

The multiple-breath washout system is designed to analyse respiratory 
gas concentrations and measure the gas volumes. Measuring the helium gas 
concentration instantly, accurately and at relatively low cost is a challenge. 
In this study, the MBHW measurement is based on two physiological as- 
sumptions: 

• Helium gas diffusing through the blood vessels is considered to be zero. 

• The gas exchanging rate(d[02]/d[C02]), i.e., oxygen consumption ver- 
sus carbon dioxide production, remains steady for different breaths and 
the existence of helium gas will not infiuence the gas exchange rate. 

The second assumption is based on our measured results of two linear 
relations between oxygen concentration and carbon dioxide concentration in 
each tidal expired breath before inhaling heliox gas. Also, the concentration 
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Figure 4: Multiple-breath-helium washout (MBHW) curve. It consists of two stages- 
washing in helium and helium washout. In helium wash-in stage, the heliox gas is washed 
into the lungs and washed out the room air. In the helium washout stage, the lungs are 
full of heliox and washed out by room air. 
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of carbon dioxide remains indenpendent of the existence of helium gas which 
has been observed from our measurements. 

The MBHW measurements can be divided into wash-in and washout 
parts. During the first part, wash-in, hehox gas (21% oxygen, 79% helium) 
is inhaled into the lungs and the air is washed out of the lungs breath by 
breath. During the second part, washout, heliox gas is washed out of lungs 
by room air. Before the wash-in, a few normal tidal breaths are taken to 
measure the gas exchanging rate. 

2.2. System Design 

The MBHW system monitors the hehum tracer gas concentration and the 
respiratory flow rate breath by breath. Thus the system can be devided into 
a gas flow measuring sub-system and the gas analysing sub-system as shown 
in figure 5. The gas flow measuring sub-system contains two flow meters to 
measure the inspired and expired gas flow with a pressure sensor to monitor 
the direction of the flow at the mouth. 

A small amount of respiratory air is pumped through the gas analysing 
sub-system, which mainly consists of a QTF gas density detector, and an 
infrared (IR) CO2 detector. Since the expired gas mixture from our lungs is 
saturated with water vapour [H2O] which interferes with the QTF and IR 
detector signal, the water vapour is removed by passing through an ice bath 
and reduced to around 0.6 % of the total gas concentration, according to the 
Goff-Gratch (1946) equation. No moisture could be detected with a commer- 
cial moisture meter (The Shaw Automatic Dewpoint Meter, Shaw, UK). The 
QTF resonance frequency (/c) gives a big response to helium gas (with about 
4 Hz frequency shift for heliox) that is related to the gas density (pg) [? ? 
]. The concentration of carbon dioxide [CO2] is monitored by an infrared 
absorption detector. A thermistor (G540 lOKl], Glass-Encapsulated Sensor, 
EPCOS, Germany) and a pressure sensor (26PGC, honeywell, Canada) have 
also been used for monitoring the gas temperature and pressure. Oxygen 
gas concentration [O2] is given by the results from the flrst few pre-washin 
breaths, and has been assumed linearly dependent on carbon dioxide concen- 
tration, based on the assumption of steady gas exchanging rate during the 
washout measurement. 

Helium gas concentration [He] can thus be given by subtracting other gas 
concentrations with two relations of 

[He] = 1 - [N2] - [H2O] - [O2] - [CO2], (10) 
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Figure 5: MBHW System. It consists of gas flow measuring sub-sytem and gas analysing 
sub-system. Gas flow measuring system consists of two flow meters and two manual vlaves 
for controlling the inhaling gas (air or heloix) . The sampling gas is taken nearby the mouth 
and then flows through a water trap which is connected with the gas analyser. The gas 
analyser mainly consists of the QTF gas density detector and a infrared CO2 detector. A 
thermistor (T) and a pressure sensor (P) are also connected to the sampling tube. 



11 



contact A contact B 



z 

L = 3.98 mm 
b = 0.60 mm 
h = 0.40 mm 
d = 0.275 mm 

Figure 6: (a) Sketch of the quartz tuning fork and (b) the main flexural resonating mode 
of 32768 Hz. 



and 

Pg = pHe[He] + PNAN2] + PHAH2O] + P0AO2] + PC0ACO2I (11) 
with [i72O]=0.61, 0< [CO2] <6, and [Osjis a linear function of [CO2]. 

2. 3. Instrumentation 

The working priciples of quartz-tuning-fork (QTF) gas density detector 
and the carbon dioxide infrared detector are described in this section. 

2.3.1. Quartz Tuning Fork Density Detector 

Piezoelectric quartz tuning forks for time-keeping element of digital watches 
are used in this study. Because of their high stability and quality factor 
(10,000 in air), quartz tuning forks have been used in various sensor applica- 
tions; for example, measuring the thermal properties of liquid helium [? ] , as 
force sensors for atomic force microscopy [??????], as SFg-gas density 
sensors [? ], and as mass sensitive sensors in liquid [? ]. The mechanical 
resonance of the tuning fork geometry is coupled electrically to contacts pat- 
terned on the surface of the quartz (figure 6) [? ]. When a tuning fork is 
placed in a fluid medium, each tine can be regarded as a damped oscillating 
cantilever [? ? ] • The equation of motion can be given by the form of har- 
monic oscillator for the fundamental mode with a driving force F at angular 
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frequency u, i.e., 

Fe^'^* = + me7^ + fcx, (12) 



d'^x dx 



where rrie is the effective mass of the cantilever, x is the relevant vibrating 
displacement, 7 is the drag force constant caused by the viscosity of the 
fluid, and k is the elastic constant of the cantilever. For a uniform cantilever 
beam with rectangular cross section, k can also be expressed in terms of the 
Young's modulus £^ by [? ] 

(13) 

where h, L, and b are the dimensions of the cantilever beam as shown in 
figure 6. The related fundamental resonance frequency / can be expressed 
in a general form of 
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(15) 
(16) 



where Cg is the geometry factor of the cantilever, pg and Pg are the density 
of the cantilever (quartz) and the surrounding medium (gas), respectively. 
For small gas density Pg compared to the density of quartz pg, the frequency 
shift (A/g) can be expressed as 

^ = ^^y^ (17) 

Jv Jv 
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where the resonance frequency in vacuum (/„) and in gas (fg) are given by 

/. = fip9-^)-^M^ (20) 
zn y Pq 



^T^^ Pq + Pg 

The resonance frequency in hehox-air mixture (///e) can be related to the 
resonance frequency in nitrogen {fN^)eind carbon dioxide (/CO2) ^ 



Afne - A/jV2 ( pHe- PN2 



fv 2 V 

A///e - AfN2 r-u ( PHe " PN2 



(22) 



(23) 



A/CO2 ~ A/iVa \PCO2 ~ PN2 , 

where pue, PN21 PCO2 ^-^^ the densities of hehox-air mixture, nitrogen, and 
carbon dioxide, respectively. For hehox (79% He + 21% O2). the frequency 
shifts about 4 Hz compared to room air. 

The net effect of the surrounding medium will cause a change in the ef- 
fective mass of the cantilever which depends on fluid density and viscosity. 
In this work, instead of measuring the frequency in vacuum {fv), the reso- 
nance frequency in pure nitrogen /jy^ and in pure carbon dioxide fc02 
measured as the reference frequencies before every the washout measurement 
to allow accurate calculation of expiratory gas density. Since the density of 
helium is much lower than of nitrogen {{pHe ~ PN2) / PN2 — —86%) compared 
to the viscosity {{pHe — Pn2)I PN2 — 5.6%), only the density effect is taken 
into account in the frequency shift. 



The quartz tuning forks (Watch Crystal NC26, Fox Electronics, USA) 
are commercially supplied in a sealed can within a small amount of helium 
gas. In this work, the quartz tuning fork is used as an actuator driven by a 
Pierce oscillator. Its resonance frequency is measured when it is exposed in 
the expired gas mixtures with the vacuum can removed. 

The QTF gas density detecting system diagram is shown in figure 7. Two 
QTFs are mounted in a plastic cell where gas passes through, one is exposed 
to the gas and the other one is left with its vacuum can in place. The one 
within its vacuum can is for compensating any change due to the temperature 
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Figure 7: QTF gas density detecting diagram. 



variation and providing a reference frequency signal. Two Pierce oscillator 
circuits have been built for each QTF and fixed on a circuit board attached 
to the QTF cell and sealed with two rubber o-rings. A commercial frequency 
detector (easyPLL plus 3.0 Detector, NanoSurf, Switzerland) with an inner 
phase lock loop (PLL) has been used for detecting the frequency shift from 
the QTF. 

2.3.2. Infrared Carbon Dioxide Analyser 

The expired gas from our breath contains around 5% of CO2. 4.2-/im 
infra-red light has been used for detecting the concentration of carbon dioxide 
which is sensitive to the CO2 but not to the water vapour [? ? ]. The infrared 
absorption by carbon dioxide follows the Beer-Lambert's law in optics which 
describes the fact that the absorption of light by a material is dependent on 
the concentraion and thickness of the material, i.e., [? ] 



where Iq is the initial light intensity while lout is the output light intensity, k is 
the absorption coefficient of the material, and / is the thickness of the material 
travelled by the light. For multiple materials, the absorption coefficient k is 
a linear combination of all compositions, i.e. 
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Figure 8: The infrared CO2 detecting cell with a temperature-control circuit. 

where Cj and are the concentrations and extinction coefficients of the ma- 
terials, respectively. 

For a small change of gas concentration, the change of the exponential 
term (e""') can be expressed in a linear form (1 — kI) and thus the output 
light intensity can be linearly related to the gas concentration. The photo- 
conductive PbSe film converts the optical energy into the electrical current 
and thus the concentration of carbon dioxide can be linearly related to the 
output voltage [? ? ]• 

As shown in figure 8, the CO2 detector contains a 4.2 fim infrared LED 
(LED42SC, Scitec,UK) and the PbSe photoconductive detector (P9696 se- 
ries, Hamamatsu, Japan) insulated from the surroundings by 1.5 cm-thickness 
of expanded polystyrene. The optical cell is temperature controlled to 49 °C 
to improve stability and to prevent the condensation of the water inside the 
cell. 

2.4- Calibration 

2.4-1 ■ QTF in Heliox-Air Mixture 

The relationship between the resonance frequency shifts and the helium 
gas concentration in air has been found by monitoring varying the known gas 
mixture composition. A linear relation is shown in figure 9 with the signal 
noise which corresponds to 0.02 Hz. 
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2.4-2. QTF and CO2 Sensor in N2-CO2 Mixtures 

During the helium washout measurements, the flowing gas mixture passes 
through the cold water trap, the QTF cell, the thermistor and pressure sen- 
sor, and the infrared CO2 cell. To get the density of the gas mixture in 
the flowing gas, a calibration has been performed by putting the 100% flow- 
ing nitrogen and 5.6% carbon dioxide gas through the gas analysing system 
before every washout measurement. Thus, the resonance frequency shifts 
A/jVj , A/CO2 ) s-iid Afg can be found. The corresponding temperature, pres- 
sure and signal from the carbon dioxide detector were also monitored at the 
same time. Figure 10 shows the singals from the QTF, infrared CO2 detector, 
pressure sensor, and thermistor with 100% N2, 5.6% CO2, and 100% CO2 
gas flowing through including values of the corresponding noise. 

2.4-3. Oxygen Concentration in Pre-washout Measurement 

Before MBHW washout measurements (before breathing in heliox) , a few 
normal tidal breaths are taken to measure the gas exchanging rate and oxygen 
concentration. The concentration of oxygen can be given by monitoring the 
gas density from the gas mitxture (N2, O2, H2O, CO2) and the carbon dioxide 
concentration with 0.61% of water vapour. Two linear phases have been 
found by plotting the CO2 concentration versus O2 concentration in each 
single expiration during pre-washout measurement shown in figure 11 with 
the corresponding signal noise. The first short phase is related to the system 
dead space where there is no gas exchange. The second linear phase has been 
defined as the gas exchanging rate in this study, i.e., the consumption of 
oxygen versus the production of carbon dioxide. The expired carbon dioxide 
concentration profile varies little from breath to breath during the washout 
measurement which implies that the existence of the helium does not affect 
the gas exchange. 

2.4.4- Signal Noise 

In this work, the infrared CO2 detector was calibrated by passing through 
different concentrations of CO2-N2 gas mixture ([CO2] from to 5.6%) with 
6mV digital noise. The oxygen concentration is given by calibrating the pre- 
washout signal with infrared sensor and QTF. The highest signal noise from 
two sensors results in 2.0% [O2] of error in room air which contains around 
21.0% of oxygen. 
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Figure 10: The output signals and the corresponding signal noise from (a) QTF fequency 
shift, (b) infrared carbon dioxide detector, (c) pressure sensor, and (d) the thermistor 
when 100% nitrogen, 5.6% carbon dioxide, and 100% carbon dioxide flowing through. 
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[CO^] (%) [CO^] (%) 

Figure 11: Systematic error caused by the sensor signal noise from two main calibrations. 

(a) is the concentration of carbon dioxide versus output infrared signal with 6mV noise. 

(b) is the oxygen concentration calibrated with both QTF and optical sensors with 2% 
[O2] error from a normal breath. 
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Figure 12: (a) Lung clearance result from a mild asthmatic female with two fitted expo- 
nential curves. The lung is occupied by 67% of a better- ventilated compartment {FRCi, 
thick curve) and 33% of less- ventilated part {FRC2, thin curve), (b) Lung clearance result 
from a normal female an exponential curve. The lung is occupied by only one compartment 
{FRC'i, thin curve) from the results of fitting with two exponential curves. 



3. Results 

3.1. Lung Clearance 

The lung clearance results from a mild asthmatic and a normal subject 
are shown in figure 12 with two fitted exponential curves representing com- 
partments with better and poorer ventilation. The result from the normal 
subject shows a nearly ideal trace with one ventilation rate. However, for 
the mild asthmatic subject, two compartments with very different ventilation 
rates that can easily be differentiated and occupy 33% and 67% of the lungs. 

The peak expiratory fiow (PEF) has been measured as a clinical refer- 
ence with a Wright metre (Mini- Wright, Standard Range, EN13826, Clement 
Clarke International, UK) on all volunteers with the predicted values from 
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the metre data shee10. The PEF measurement has been regarded as a sen- 
sitive technique for identifying airflow obstruction [? ]. From our MBHW 
results, the measured lower PEF values are not necessarily related to the 
lung obstruction but somehow reflected in a larger proportion of poorly ven- 
tilation compartment {Vpoor/ PRC) with very different ventilation rates. For 
example, the PEF values for two male subjects (No.l and 3) are much smaller 
than the predicted values, and both of them have very different ventilation 
rates. However, the measured PEF value for the asthmatic subject (No. 11) 
is higher than the predicted value which can be a result of regular exercise 
habbits, while both of the compartments and ventilation rates are apparently 
different. 

3.2. Normalised Phase III Slopes 

Figure 13 shows the results from two mild asthmatics and two normal 
subjects in which the normalised phase III slopes of helium during (a) helium 
washin and (b) helium washout process are plotted versus the turnover. In 
(a), the wash- in results from a mild asthmatic female shows both of Scond 
and Sacin values are higher than the normal ones which indicates the higher 
ventilation inhomogeneities in both conductive and acinar airways. In (b), 
the washout results from a mild asthmatic male shows a higher Scond value 
than the normal one but similar Sadn values which indicates the asthmatic 
male has no apparent higher ventilation inhomogeneity in the acinar airways 
than the normal volunteer. 

The Scond and Sadn values from the all the volunteers have been listed 
in table 2. The results from four groups (normal, normal with asthma and 
COPD history, smoker, mild asthmatic) are compared and plotted in figure 
14 and 15. The mild asthmatic group shows higher Scond and Sadn values in 
wash-in but not much difference in washout measurements compared to the 
rest of the groups. For mild asthmatics, subject No. 6 has higher Scond values 
while subject No. 7 has higher Sadn values, and No. 11 has higher values in 
both Scond and Sadn compared to the normal group. For the two smokers, the 
Sadn values are higher compared to the normal-subject group which reflects 
with ventilation inhomogeneity in the acinar airways. 
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Table 1: MBHW Results - Lung Clearance 



Gender, 
Age(yrs) 



Height Weight 
(cm) (kg) 



PEF'^/PEFP 
(%) 



Qr FRC Vr 
(L) (L) 



^better / FRC Vpoar / FRC 



Rbetter/R'i 



poor 



LCI He 



M, 25 
M, 24 
M, 25 
M„,„ 27 
M,h, 29 
M„, 29 
M„, 47 
M„ 27 
F, 26 
F, 26 



Fa, 24 



F„ 35 



173 
178 
195 
175 
185 
165 
175 
175 
166 
155 
155 
163 



100 
66 
75 
65 
78 
59 
76 
74 
50 
57 
54 
53 



79.77 

95.04 
66.72 

106.80 
91.81 
75.63 
91.20 
129.03 
95.24 
99.30 
107.14 



1.50±0.13 
1.65±0.06 
1.30±0.11 
1.64±0.09 
2.03±0.08 
1.97±0.11 
1.74±0.09 
1.52±0.07 
2.33±0.08 
1.90±0.11 
2.39±0.08 
1.60±0.07 



4.16 
2.76 
5.77 
3.61 
3.58 
2.64 
3.71 
3.40 
3.01 
1.87 
2.03 
1.86 



1.02 
0.91 
2.09 
0.92 
1.49 
0.56 
1.52 
0.88 
0.89 
0.69 
0.49 
0.73 



0.41 
0.93 
0.24 
1 

0.84 
0.49 
0.68 
0.99 

1 

0.98 
0.33 

0.84 



0.59 
0.07 
0.75 

0.16 
0.51 

0.32 
0.01 

0.02 
0.67 

0.16 



I. 8421 
2.9174 
2.2336 

1 

1.6648 

II. 0674 

1.5007 
1.5470 
1 

1.4131 
4.6002 

1.6829 



5.64 
5.15 
7.70 
4.09 
6.30 
6.95 
6.69 
5.04 
4.67 
5.93 
5.88 
6.85 



Lung Clearance Results from four male (one with asthma history M^h, No. 4*, and one with COPD history, 
Mch, No.5*) and two female normal volunteers, two male mild asthmatics (M^, No. 6"^ and 7+) and a 
female mild asthmatics (F^, No. 11+), and a female (F^, No.l2**) and male (M^, No.8**) smoker. Two 
normal male subjects (No.l and 3) having smaller ratio of the measured PEF values (PEF™') over the 
predicted values (PEF^) are found to have larger capacities of poorly ventilated volumes Qr is the 

gas exchanging rate defined in section 2.1. 




Figure 13: (a) The normalised phase III slopes of helium during helium wash-in process 
from a mild asthmatic and a normal female volunteer. The results from the asthmatics 
show higher absolute values of Scond (S^ond =0.016±0.004) and Sactn (S^ctn = 0-41 ± 
0.0007) compared to the normal ones (S';^„„rf=0.0055±0.0015, \S2c^n\ = 0.12 ± 0.0003). 
The arrow indicates the normalised phase III slope of the first breath from the asthmatics 
is much steeper compared to the normal volunteer which results in a higher absolute 
value of Sacin and implies higher ventilation inhomogeneity in the acinar airway, (b) 
The normalised phase III slopes of nitrogen during helium washout process from a mild 
asthmatic and a normal male. The results from the mild asthmatic show a higher Scond 
value but similar Sacin value compared to the normal ones. 
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MBHW results from all volunteers 



No. 


Gender, 


LCI 


qHe,in 

acin 


qHe,out 

acin 


QHe,in 

COTid 


oHe^out 

cond 




Age(yrs) 




*io-3(L-i) 


*10-3(L-i) 


*10-3(L-i) 


*10-3(L-i) 


1 


M, 25 


5.65 


-154.35±0.71 


131.61±10.70 


2.20±0.85 


70.98±13.8 


2 


M, 24 


5.15 


-195.30±0.22 


46.42±4.20 


5.75±0.52 


56.04±7.58 


3 


M, 25 


7.70 


-217.30±0.65 


33.80±1.60 


6.64±1.66 


23.09±2.57 


4* 


Mah, 27 


4.09 


-277.79±2.00 


139±11.10 


9.75±5.95 


40.57±20.91 


5* 


Me,., 29 


6.29 


-115.40±2.24 


82.30±2.50 


2.75±0.39 


16.24±2.87 


6+ 


M„, 29 


6.95 


-196.43±0.61 


135.00±0.98 


7.39±1.91 


95.91±5.21 


7+ 


Ma, 47 


6.68 


-240.99±0.53 


239.10±2.10 


6.01±0.99 


26.79±3.10 




M„ 27 


5.04 


-242.20±0.24 


181.40±5.70 


3.88±0.75 


45.89±12.68 


9 


F, 26 


4.67 


-120.10±0.38 


24.51±2.09 


5.49±0.75 


74.04±26.65 


10 


F, 26 


5.93 


-118.80±0.51 


128.10±36.10 


3.59±0.55 


75.54±22.52 


11+ 


F„, 24 


5.88 


-401.60±0.78 


228.30±25.30 


16.45±2.05 


109.16±43.52 


12** 


F«, 35 


6.85 


-377.40±0.95 


376.30±12.3 


14.06±2.59 


30.29±21.55 



Table 2: MBHW Results from all volunteers. Scond values during helium wash-in 
and washouts. The lung clearance index from helium washout results are also 
listed. 
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normal 



(b) 



asthma and 
COPD history 



normal 



asthma and 
COPD history 



smoker 



smoker 



asthma 




asthma 



Figure 14: MBHW results of Scond values from all volunteers into four groups where the 
errors have taken the standard deviation for each group with the mean values, (a) is from 
the helium wash-in and (b) is from the helium washout results. Mild asthmatic group has 
shown a slightly higher mean value than normal group. 



3. 3. Smokers 

An extended phase II and delayed phase III have been observed in most 
of the single breaths from four smokers as shown in figure 16. For smokers 
with longer smoking history (>10 years), this feature becomes more appar- 
ent which implies higher uneven ventilation between conductive and acinar 
airways. The curved transition to phase III, which is little documented in 
nitrogen washout measurements, results in a difficulty of defining phase III 
but may imply an early stage of phase III steepening which has been ob- 
served in COPD patients. Also, small and unsteady tidal volumes have also 
been observed from four smokers compared to the normal volunteers with 
similar FRC The helium washout results from two of the smokers have to 
be excluded because of their constant coughing or short of breath during the 
measurements (10- and 20- year smoking history). The washout results from 
the other two smokers (5- and 15-year) have included the curved transition 
into phase III. 



This can be found on the website http://www.peakflow.com/top-nav/norrnaLvalues/index.html 
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(b) 0.4 
0.3 

I_i 

0.2 

o 

w 

0.1 




normal 



• He washout 



normal 



asthma & 
COPD history 



smoker 



asthmatic 



asthma & 
COPD history 



smoker 



asthmatic 



Figure 15: MBHW results of Sadn values from all volunteers into four groups where the 
errors have taken the standard deviation for each group with the mean values, (a) is from 
the helium wash-in results and (b) is from the helium washout results. Results show that 
smokers have higher Sacm values in both wash-in and washouts. 
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Figure 16: The single breath curves from MBHW measurements with four smoker volun- 
teers having different length period of smoking history compared to the normal subjects. 
For smokers there is a more extended phase II which causes the phase III linear region less 
clearly defined. For smokers with longer period of smoking history (>10 years), the feature 
becomes more evident which implies a higher uneven ventilation between conductive and 
acinar airways. 



28 



4. Conclusion 

The MBHW system containing a QTF gas density sensor and a infrared 
CO2 detector has been developed for MBW lung function study. The system 
gives accurate values of gas concentrations and can be used to extract the 
LCI as well as Scond and S'^cm-The MBHW results from 12 volunteers show 
a difference between normal and asthmatic (or smokers) subjects. Subjects 
with lower PEF values (compared to the predicted values) are found to have 
smaller portion of better- ventilated lung capacity, which is not necessarily 
related to lung disease. 

Higher Scond values or Sadn values have been found on mild asthmatics 
compared to normal subjects which reflect the inhomogeneis ventilations in 
mildly diseased lungs. The higher Sadn values have also been found on the 
smoker with 15-year smoking history. 

The results from the wash-in process show a more apparent difference in 
Scond and Sadn between the four groups indicating that it may be a more 
sensitive manuevour compared to the washout process. 

An extended phase 11 and delayed phase III have been observed in single 
breath curves from four smokers with different length of smoking history. 
For smokers with longer smoking history (>10 years), this feature becomes 
more apparent which implies higher uneven ventilation between conductive 
and acinar airways and may indicate the early stage of phase III steepening 
observed from COPD patients. 

Appendix A. Lung Clearance 

The washout lung clearance curve is taken by plotting the mean con- 
centration of gases (helium or nitrogen) as a function of turnover TO. The 
mean concentration of gases (c„ from nth breath) have been taken as the 
expired amount of gas (/q^^ CndV) divided by the expired tidal volume (Vr) 
from each breath. The turnover TO is the accumulative expired tidal volume 
(S ^t) divided by the functional residual capacity (FRO). In the helium or 
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nitrogen washouts, FRC has been simply derived from 



FRC - Co ... 
^ FRC + Vry ^^-'^ 
FRC ■ ci , , ^, 



FRC + Vt, 



(A.3) 
(A.4) 



Co ~ c„ 

where c„=o,i,2,---n-i and Vr,„(„=o,i,2,---n-i) a-i'e the mean concentration of gas 
and tidal volumes from the initial gas in the lungs, the first, second, and 
{n — l)th expiration, respectively. This equation is based on the assumption 
of well-mixed gases in the lungs. Since the FRC has been defined as the 
lung volume at the end of expiration, it is dependent on the expired tidal 
volumes and is slightly different breath by breath. In this study, only gas 
concentrations higher than 0.1% have been taken into acount for the FRC 
calculation. 

On the basis of the two-compartment model, human lungs can be ab- 
stractly regarded as two parallel compartments with different ventilations [? 
]. In this section, the washout results from the helium washout process will 
be discussed. The mean concentration of helium from each breath is plotted 
as a dual-exponential function (the summation of two exponential functions) 
of turnover (cumulative expired tidal volumes divided by FRC). It can be 
expressed as 

FRC = FRC1 + FRC2, (A.6) 
Vt = V1 + V2, (A.8) 

+ = frc.Tfrc, ''-^ frc^Tfrc, '''- ^^-'^ 



Cl,n — CaO 



TO 



TO 



^ '^'{ frc!7vJ ' ^^-''^ 
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where FRCi and FRC2 are the residual capacities of two compartments 1 
and 2, Ci „ and C2,n are the mean hehum concentration from the nth breath, Vi 
and V2 are the expired tidal volumes of two compartments. If compartment 1 
represents the better- ventilated compartment compared to the compartment 
2, then 1 has a higher washout rate, e.g., j.^^'^Vi frci+Vi ■ ratio of 
represents the proportion of the better- ventilated part of the lung. 

In this work, breaths from the first expiration to the end (> 0.1%) of 
helium gas concentration during helium washout have been taken into ac- 
count for this two-exponential curve fittings. The intial gas concentraion 
C(jo is thus the mean gas concentration from the first breath. If the initial 
gas concentration was taken as the gas concentration in the lungs before the 
washout measurement, the lung clearance washout curve was unable to be 
fitted well with a sum of two exponential curves. This is likely to result from 
the imperfect gas mixing such that the mean concentration of the expired gas 
is smaller than the mean concentration of gas from the whole lungs. FRCi, 
FRC2 are given by Eq. A.9 and Vi, V2 are given by Eq. A.IO and A.ll. 

For more complex models with 50 compartments can be found in the 
studies by Lewis et al. [? ? ]. A model of tidal emptying pattern has 
also been developed in a few modern studies by Whiteley et al. with the 
consideration of gas re-breathing from the dead space [???]. The math- 
ematical model based on the one- and 50-compartment model are described 
and compared with the results from the inert gas elimination technique. Un- 
like the MBW technique, the tracer gas is inserted into the blood and the 
concentration of the tracer gas is sampled by taking the blood sample. From 
the inert gas elimination results, the breathing pattern is more likely tidal 
emptying instead of continuous which means the gas re-breathing from the 
dead space is unable to be ignored. 
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